Fano-like interference of plasmon resonances at a 
single rod-shaped nanoantenna 



F Lopez- Tejeira, R Paniagua-Dommguez, 
R Rodriguez-Oliveros and J A Sanchez-Gil 

Institute de Estructura de la Materia (lEM-CSIC), Conscjo Superior de 
Investigaciones Cientfficas, Serrano 121, E-28006 Madrid, Spain 

E-mail: flt@ieiii.cfmac.csic.es 

Abstract. Single metallic nanorods acting as half- wave antennas in the optical range 
exhibit an asymmetric, multi-resonant scattering spectrum that strongly depends on 

both their length and dielectric properties. Here we show that such spectral features 
can be easily understood in terms of Fano-like interference between adjacent plasmon 
resonances. On the basis of analytical and numerical results for different geometries, we 
demonstrate that Fano resonances may appear for such single-particle nanoantennas 
provided that interacting resonances overlap in both spatial and frequency domains. 
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1. Introduction 

Different experimental PP-[S] and purely theoretical [S]-[in] investigations have shown 
that metallic nanorods act as standing-wave resonators for locahzed plasmon resonances 
in the optical regime, thus exhibiting geometrical half-wavelength resonances with 
spectral positions depending mainly on the length of the rods. This particular type of 
so-called "optical nanoantennas" have raised the prospect of significant improvements in 
fields such as photodetection [TT], field-enhanced spectroscopy [12] or control of emission 
direction in single-molecule light sources [T3] . 

Generally speaking, most of device-oriented studies are focused on nanoantennas 
operating at the dipole-like resonance. However, structures with a high aspect ratio 
may support additional resonances that have usually been the subject of a more 
fundamental research work. Hence, several authors have already elucidated the scaling 
properties of high-order longitudinal modes, as well as their dependence on shape, size, 
orientation and dielectric environment by means of diverse approaches and techniques 
[21 E], [I1]-[IZ]- Nevertheless, a relevant issue has yet to be addressed for multi- 
resonant nanoantennas, that is the emergence of asymmetric line profiles in single- 
particle extinction or scattering spectra. Interestingly, such a feature seems to go almost 
unnoticed for the nanoplasmonics community, despite being apparent in some previous 
references [21 [TJ [ISl [17] . In fact, to the best of our knowledge, the only explicit report on 
the occurrence of Fano-like asymmetric line shapes in the scattering spectra of a single 
silver nanorod can be found in a recent paper by Reed et al. [18|, though the emphasis 
is put on dimer structures therein. 

In this work, we show that these asymmetric line profiles can be easily understood in 
terms of the so-called Fano-like interference between localized plasmon resonances that 
has been recently reported for a variety of coupled metal nanoparticles [l9]-[22]. Being 
more precise, we present a simplified analytical model that describes spectral features 
of a single rod-shaped nanoantenna in terms of Fano-like interference. Contrary to 
the common assumption that interference does not play any role in total scattering or 
extinction of a single metallic surface, we find a good agreement with numerical results, 
which are attained through the separation of variables (SVM), finite element (FEM) and 



surface integral equation (SIEM) methods (see Appendix A for a succinct description 
of calculation techniques). Furthermore, we make use of explicit expressions for light 
scattering by spheroids to conclude that not only spectral but also spatial overlap (i.e. 
non-orthogonality) between interacting modes underlies the emergence of such single- 
rod resonances. This points out the need of being extremely cautious when applying the 
premises of standard Mie theory to particles that significantly depart from sphericity. 
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2. Single metallic nanorods acting as half- wave nanoantennas 
2.1. Fano-like interference of longitudinal plasmon resonances 

Let us begin by briefly reviewing tlie basics of liglit scattering by a single metallic 
nanorod with dielectric function Em that is surrounded by a medium with constant 
permittivity Ed'- On the assumption that rod diameter D is much smaller than its total 
length L, the electromagnetic response to p-polarized light impinging perpendicular to 
the long side is fully governed by longitudinal modes. The fundamental resonance Arel 
can thus be described as a dipolar excitation of charges at the rod surface, with its 
wavelength exhibiting a linear dependence on L , Aril oc 2L. For a perfectly conducting 
material of negligible thickness, A^-el is precisely equal to 2L, whereas Alel ^ 2L at 
optical frequencies [11 El E] • As rod length increases, additional resonances may appear. 
Following Khlebtsov and Khlebtsov [7], we assume the position of any longitudinal 
resonance to be described by the following approximate scaling law 

^r:>B, + B,^ (1) 

where n is an odd integer and Bq, Bi are coeflicients that depend on the actual 
geometry and dielectric environment of the system and have to be determined from 
linear regression. With respect to Bq, Bi, we also have to notice that, although not used 
in this paper, explicit expressions can be obtained within the framework of Novotny's 
model for effective wavelength scaling at optical antennas [8|. 

As mentioned in the previous section, the interaction of adjacent resonances has 
been found to be compatible with a Fano-like interference model |23],[20]-[25], where 
the lower resonance plays the role of continuum in canonical Fano line shape [26] . Given 
that light scattering by a rod of flnite length cannot be treated in a closed form within 
the framework of standard Mie theory [2^ , we assume a heuristic line shape for the first 
approach. Hence, 

Qsca^l/Mr (2) 

with 



f{u;)=A{u;) + F 



{uj - UJi) + ihi {uj - Ws) + ih-i 



(3) 



where A{(jj) is a slowly varying amplitude and F stands for the complex amplitude of 
the fundamental resonance, which is described in terms of its real central frequency Ui 
and spectral width hi. Dimensionless real parameter q modulates the interaction with 
the adjacent plasmon resonance (analogously defined by and 63), thus governing the 
line shape asymmetry. 



2.2. Prolate spheroidal nanorods 

From the computational point of view, the shape of half-wave nanoantennas is 
usually modeled by a right circular cylinder with either fiat or hemispherical ends 
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Figure 1. (a) Calculated scattering efficiency as a function of wavelength for a single 
Ag spheroid surrounded by glass (e^ = 2.25). Incident field is p-polarized and impinges 
perpendicular to the rotation axis of the spheroid. Different curves correspond to 
increasing values of L, whereas D is set to 30 nm for all calculations, (b) Linear 
scaling of resonant wavelengths in panel (a) as a function of normalized aspect ratio 
L/nD. Dashed line marks the best fit to Eq. Q. 
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Figure 2. Calculated scattering efficiency (solid line) as a function of photon energy 
for a single Ag spheroid (L = 340 nm; D = 30 nm) in — 2.25. Incident field is 
p-polarized and impinges perpendicular to the rotation axis of the spheroid. Dashed 
line renders the best-fitting curve to Eq. ([2| . Obtained values of bj\ , 0J3 and q are also 
shown. Lower right inset panel depicts the obtained values of q for every L in figure [T] 



(spherocylinder). However, we have found the prolate spheroid to be the most convenient 
geometry to start with, because of the following reasons: (i) Previous works [2H1 [22] on 
the basis of SVM have provided us with a very efficient approach to calculate extinction, 
scattering and absorption cross-sections even for very elongated nanoantennas; (ii) In 
addition to its low numerical cost, such a theoretical framework also makes apparent 
the origin of unexpected asymmetry in line profiles, as detailed hereafter; (iii) Prolate 
spheroids are not only of academic interest, as they accurately describe the so-called 
nanorice structures [T71 130] . 

In figure [T] we present the calculated scattering efficiency Qsca for a single silver 
spheroid surrounded by glass {sd = 2.25) under the assumption that incident field 
is |}-polarized and impinges perpendicular to the long side of the rod. Different curves 
correspond to increasing values of total length L within the [100, 400] nm range, whereas 
the polar diameter D is set to 30 nm for all calculations. As can be seen, the position 
of resonances increases linearly with L and such displacement is fairly well described by 
([1]) (see figure [T]^b)). For L/D > 5, the peaks arising from resonances with n = 1 and 
n = 3 are clearly apparent, as it is the asymmetry of the line shape between them. 

In order to test out our approach, let us take a closer look to the curve corresponding 
to 340 nm-long spheroid. Figure [2] renders the calculated scattering efficiency as a 
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function of energy and its best-fitting curve to (|2]). It may be seen that our heuristic 
Fano-hke hne shape agrees very well with the full electromagnetic calculation. Besides, 
the obtained value q = 0.82324 is consistent with asymmetric profiles being described 
by |g| ~ 1 [19]. Such an agreement extends to the whole range of rod lengths in figure [l| 
as summarized in the inset panel (See Appendix B). 

For the case of oblique incidence (i.e. if the angle a between the incident light vector 
kinc and the rotation axis of the spheroid is not equal to 90°), "even" modes become 
accessible [21 |7l [I7| and should therefore be incorporated to our analysis. Different 
curves in [sj^a) correspond to the calculated Qsca for decreasing values of a and the same 
geometry, polarization and dielectric environment as in figure |2} Logarithmic scale is 
used for optimal visualization of less intense features. As can be seen, the oblique line 
shapes differ quite little from that for normal incidence at the vicinity of resonances 
with labels n = 1,3, except from that the intensities of "odd" peaks exhibit a linear 
dependence with sin^ a (see lower left inset panel). However, an extra resonance develops 
just on top of base line at u; ~ 1.045 eV. Given that it is located within the [ci;i,Ci;3] 
range, we label it as = 2. This new peak reaches its maximum for a = 45° and seems 
to be symmetrical with respect to its central frequency. Hence, we can conjecture that 
symmetry precludes interference between longitudinal modes with different parity, so 
that only an additive term accounts for the contribution of this resonance to Qsca, 

Qscaioo, a ^ 90°) ^ |/(a;)r + 'f' (4) 

where F2,uj2, &2 are the complex amplitude, the central frequency, and the spectral width 
of the resonance with label n = 2, respectively. 

Comparison between calculated Qscai^^, ce = 45°) and its best-fitting curve to Q in 
figure |3](b) confirms our hypothesis of the lack of interaction between "odd" and "even" 
resonances. As summarized in the upper right inset panel, the only effect of oblique 
incidence (aside from the emergence of a new peak!) is the already mentioned global 
quenching of Qsca{(-^res^^^) (i.e., the intensities of "odd" peaks) and therefore that of q. 
Being the understanding of line shape asymmetry the main goal of our work, we do 
not discuss on oblique incidence any further, except for the detailed description of the 
dependence of different parameters on a that is presented in [Appendix B However, 
we have to point out that inter-parity coupling may be allowed for a symmetry-broken 
configuration, such as depositing the nanospheroid onto a dielectric substrate, as recently 
proposed for the nanocube geometry [31]. In fact, we can easily envisage that the two 
different intra- and inter-parity mechanisms' operating simultaneously opens a very 
interesting scenario that certainly warrants further investigation. 

Going beyond heuristic description brings us up against the actual meaning of 
From a formal point of view, it is clear that f{uj) plays the role of effective polarizability 
in a somehow generalized Rayleigh-Gans theory for nonspherical particles. Nevertheless, 
there is still the concern of how to explain the emergence of asymmetry in line profiles. 
As previously mentioned, Fano resonances require an observable that is sensitive to 
interference, but standard Mie theory predicts the total scattering by a single metallic 
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Figure 3. (a) Calculated Qsca as a function of photon energy for a single Ag spheroid 
(L = 340 nm; 13 = 30 nm) surrounded by Ed = 2.25. Incident field is p-polarized and 
impinges with angle a with respect to the rotation axis. Different curves correspond to 
decreasing values of a. The intensities of "odd" peaks as a function of a are presented 
at lower left inset panel, (b) Comparison between calculated Qsca and its best-fitting 
curve to Eq. Q for a = 45°. Obtained values of u}i,uj2,i^3 and q are also shown. 
Upper right inset panel depicts the obtained values of q for every a in (a). 
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surface to be proportional to the mere sum of intensities. 

Such a discrepancy can be easily explained for the case of prolate spheroidal 
particles by means of the SVM formalism. According to Reference |29], the scattering 
efficiency of a prolate spheroid for p-polarized light impinging at normal incidence is 
given by 

. f OO OO OO OO p 



1=1 1=1 m=l n=m 



+ ( 6L1 [a^iy^t\c.. Q) - (6Lt)*47Hc„ c.) ) + (5) 



Without entering into further details, let us say that t/™^ and wj'™^ are different 
integrals of the normalized prolate angular spheroidal wave functions with normalization 
coefficients Nij{cd), where q = ^1?- — D"^, = 27r/A and A is the wavelength 



of incoming (or outgoing) radiation. More importantly, a-^^ and b\'^\ hf^ represent the 
expansion coefficients of the azimuthal components of scattered electric and magnetic 
fields respectively. Single or double subscript marks whether the expansion accounts 
for the axisymmetric (i.e. the one that does not depend on the azimuthal angle) or 
non-axisymmetric (i.e. azimuth-dependent) part of the corresponding field. For the 
axisymmetric part, the introduction of Debye's potential representation of vector fields 
(i.e. the potentials that solve the light scattering problem for spheres) leads to the 
term following the curly bracket in ([s]), which closely resembles Mie's result. However, a 
proper representation of the non-axisymmetric part of the electromagnetic fields requires 
a combination of Debye's potential with that used to deal with light scattering by an 
infinitely long cylinder, namely Hertz's [32]. Hence double-subscript interference terms 
arise in Qsca-, thus providing a mathematical description for the interaction between 
adjacent resonances in prolate spheroids. Note that such interference terms are not 
negligible only if the interacting fields are non-orthogonal, that is, spatial overlap 
between interfering plasmon modes is required for Fano-like interference. 

Although Eq. ^ is only valid for spheroids, we conjecture the underlying sphere- 
like vs. cylinder-like interference mechanism to operate for any single nanoantenna with 
a similar geometry. In order to determine if there is any real substance in our guess, 
we present in figure |4] the calculated scattering efficiency for a single 550x30 nm Ag 
nanoantenna embedded in Sd = 2.25 assuming three different geometries: a circular 
cylinder with flat ends, a spherocylinder and a prolate spheroid. Numerical values are 
obtained from either FEM or SVM. Logarithmic scale is used for the sake of a better 
comparison. As can be seen, the three curves can be fltted to Eq. ^ even for the case 
of a flat-ended cylinder, although only a modest g ~ 0.15 is obtained for such a high 
aspect ratio. Hemispherical ends result in more than a 30 percent increase of asymmetry 
parameter, which rises to its maximum value of g ~ 0.92 for prolate spheroidal geometry, 
where spherical and cylindrical features interact in the most efficient way. These results 
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Figure 4. Calculated scattering efficiency as a function of photon energy for a single 
Ag nanorod {L = 550 nm; D ~ 30 nm) surrounded by glass (e^ — 2.25). Incident 
field is p-polarized and impinges perpendicular to the rotation axis. Black, red, and 
blue curves correspond to flat-ended cylinder, spherocylinder and prolate spheroid, 
respectively. Open symbols render the best-fitting curves. Obtained values of q are 
also shown. 



are consistent with our previous formal discussion on potential theory and point out the 
subtle balance of different contributions to light scattering. 



2.3. Spherocylinder- shaped nanorods 

In figure |5] we present a separate plot of scattering efficiency as a function of photon 
energy for the spherocylinder-shaped nanorod in figure |4j As already mentioned, the 
asymmetric line profile can be fitted to (g (dashed line). The obtained q = 0.20294 
accounts for the moderate interaction between adjacent resonances located at Wres = 
0.45563 eV and Urll = 1.19918 eV. Please notice that spectral features at u > 1.7 
eV suggest the need to include subsequent resonances in ([s]), although that refinement 
is beyond the scope of our present work. Inset panels at the right hand side show the 
calculated Qsca curves (upper) and their corresponding q values (lower) for different rods 
with D = 30 nm and L within the [340, 550] nm range. As expected, the interaction 
between adjacent resonances is significantly lower than that of prolate spheroids with 



the same L, D parameters. (Details on fitting are presented in Appendix B ) 

The normalized electric near-field distribution in the xz plane at most significant 
values of photon energy is presented in the three panels at the left hand side of figure [5] 
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Leaving aside the typical ra-node quasi-standing-wave patterns at oorei and ujres, let us 
concentrate on the plot for u = 1.12139 eV. Given that Qsca reaches its local minimum 
within at this precise value, we expect some destructive interference to appear in spatial 
domain. Such interference pattern can be noticed for the field distribution inside the 
volume of the rod, which is limited by dashed lines. For the sake of clarity, we present in 
the bottom panel the line profiles at x = (dotted lines on contour plots). In short, for 
u = Umin, interference cancels a significant part of the field intensity that is present at 
the zones marked by descending arrows for lo = oorei- On the other hand, field intensity 
at the central part of the rod (ascending arrow) for u = Umin is slightly enhanced with 
respect to that for uj = Ur^l, thus resembling the field pattern corresponding to 



2.4- The simplest case: a rectangular nanowire 

For a better understanding of the elusive spatial overlap between modes, we now consider 
light scattering by an infinitely long rectangular nanowire with large aspect ratio, which 
seems to be the simplest geometry that clearly shows such behavior. In figure |6] we 
plot the scattering efficiency (blue curve) and the near-field amplitude (black curve) at 
normal incidence calculated using SIEM for a silver rectangular nanowire with L = 600 
nm and D = 10 nm surrounded by = 1. Near- field amplitude is now evaluated 
3 nm outside the end of the nanowire and normalized to the incident field at this 
point. However, the two curves are then re-normalized to unity to more clearly show 
their spectral shift. Such a change in the position of maxima, which has recently been 
explained in terms of driven and damped harmonic oscillators |33], will be relevant for 
our discussion. 

In this simple geometry, the very broad dipole-like resonance in Qsca at A ~ 1500 
nm is more similar to the "canonical" Fano continuum than those in the previous 
configurations and therefore provides a strong spectral overlap with the narrower 
3A/2-mode. With respect to spatial overlap between different longitudinal plasmon 
resonances, we find that it can be monitored by means of the real part of the normal 
component of electric field (i.e. Enear figure [6]), which is directly proportional 
to surface charge distribution. Given that, in the electrostatic limit, surface charge 
distribution can be approximated by a sum of sinusoidal functions with argument an 
integer multiple of irx/L, we finally arrive to the following expression 

N p -| 

Re[Er7^r](^)^ Yl ^"^i^ ^ (6) 

n=2fe+l L 

for which the set of coefficients {An} can be straightforwardly determined from Fourier 
analysis. 

According to this approximation, we expect surface charge distributions at the 
maxima of near- field amplitude (and not of scattering efficiency!) to have a net 
dipolar moment and therefore to exhibit a pattern with an odd number of nodes and 
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Figure 5. Calculated scattering efficiency as a function of photon energy for a single 
silver spherocylinder (L = 550 nm; D = 30 nm) surrounded by ej, = 2.25. Incident 
field is p-polarized and impinges perpendicular to the long side of the rod. Dashed 
line renders the best-fitting curve to [2j Inset panels at the right hand side show the 
calculated Qsca curves (upper) and their corresponding q values (lower) for different 
rods with D — 30 nm and L within the [340, 550] nm range. The normalized electric 
near-field distribution in the xz plane at most significant values of photon energy is 
presented in the three panels at the left hand side. Line profiles at x = are shown in 
the bottom panel. 
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Figure 6. Calculated scattering efficiency (blue) and near-field amplitude (black) 
as a function of wavelength for a single Ag nanowirc {L = 600 nm, D = 10 nm) 
surrounded by e^; = f. Incident field is p-polarized and impinges perpendicular to 
the xy plane. The spectra have been normalized to unity. Inset panels show the 
calculated Re[£'"^'"'/i?™'^] along the nanowire (solid) as well as its fitting to [g] for 
X S [—L/2,+L/2] (dash-dotted), evaluated at the three most significant wavelengths 
in spectrum. Vertical dashed lines mark the limits of the nanowire. 



charges of opposite sign at its ends. Surface charge distribution will gradually modify 
as wavelength varies and adjacent resonances come into play. In fact, somewhere 
within the intermediate region, destructive interference between modes reduces near- 
field amplitude to a minimum, as can be seen in figure |6| In order to quantify such 
features we find it useful to define the following auxiliary magnitudes: as a measure of 
the contribution of n-resonance at a given wavelength, we introduce its relative weight 
as 



which will vary from 1 in the vicinity of n-maximum to /„ ~ as the charge 

reaches its subsequent resonance. Additionally, we take the right hand member of ^ 
as the measure of surface charge at x = +L/2: 

N 

A„= (-i)^""''/'^n ~ Re[i?r7£^r]x=+L/2 (8) 



n=2k+l 
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Table 1. Obtained values of /i, /a, /s, A^- for inset panels in figure [6] 
A = 1993 nm A = 939 nm A = 700 nm 



/i 0.8304 0.7289 0.0785 

/s 0.08887 0.2196 0.8279 

/s 0.08075 0.05148 0.09362 

3.8737 0.01947 1.95 



Inset panels in figure [6] show the calculated 'Re[E^^"'^ / El"''^] along the nanowire 
(solid) as well as its Fourier series up to = 11 for x G [—L/2, +L/2] (dash-dotted), 
evaluated at the three most significant wavelengths in spectrum: At A = 1993 nm, 
surface charge exhibits a typical dipole-like distribution that is in good agreement with 
the obtained /i = 0.8304, /g = 0.08887, whereas /s = En>3 fn = 0.08075 accounts 
for the rest of contributions. A similar behavior can be found at A = 700 nm, but 
longitudinal mode with n = 3 now plays the main role and the values of /i,/3 are 
almost interchanged, namely 0.0785, 0.8279. For the relative minimum of near-field 
amplitude at A = 939 nm, the situation is completely different. Although obtained 
values of /i, f^, /s are not too different of those at A = 1993 nm (see table [T]), the field 
profile resembles that of A = 700 nm. However, surface charge reduces drastically at 
the ends of the nanowire, which is consistent with the expected minimum in dipolar 
moment. In fact, the obtained A^- = 0.01947 is two orders of magnitude less than those 
for A = 1993, 700 nm (see last row in [T| and nearly approaches to the condition of 
complete destructive interference at x = ±L/2 in (|6]), which is given by A^- = 0. We 
find this direct evidence of spatial overlap between interfering plasmon modes to confirm 
the consistency of our Fano-like model when applied to a generic nanorod geometry. 

3. Conclusions 

In conclusion, the emergence of asymmetric line profiles in the total scattering spectra 
of single metallic nanorods acting as half-wave nanoantennas is explained in terms of 
Fano-like interference between adjacent, odd-order plasmon resonances. Such a feature 
can be understood as originated from the interplay between the different contributions 
that properly describe the non-axisymmetric part of the scattered field. Our analytical 
and numerical results show that Fano resonances can be excited on diverse geometries 
(elongated nanospheroids, nanorods and nanowires), provided that interacting modes 
overlap not only in energy but also spatially (non-orthogonality). This finding makes 
single-particle nanoantennas especially suitable for a wealth of applications where the 
sharp, environment-sensitive Fano-like spectral profile can be crucial [34j- Moreover, the 
underlying mechanism of the resulting Fano resonances has far reaching implications in 
exploring other single nanoparticle configurations designed for specific applications. 
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Appendix A. Calculation techniques 

Calculated scattering efficiencies for spheroids were attained by means of a modified 
version of the F77-code made publicly available at the web site of the Jena-St. Petersburg 
Database of Optical Constants (JPDOC) [35j. Detailed description of SVM formalism 
and benchmark results for the code can be found in References |29] and [32], respectively. 

Results for cylindrically-shaped rods with either flat or hemispherical ends were 
obtained with the RF module of COMSOL Multiphysics flnite element software, version 
4.2. The computational domain consisted on the nanorod (which, for the case of the 
spherocylinder, was modeled as the union of a cylinder, with length L — D and diameter 
Z?, and two spheres of radii R = D/2 centered in each of the cylinder flat faces) and 
two concentric spheres of radii -Ri = 0.7(L — D) and R2 = 1.1{L — D) surrounding it. 
The space between the smallest sphere and the rod models the medium in which the 
nanoantenna is embedded (glass). The subdomain comprised between the two spheres 
is set to be a Perfectly Matched Layer (PML), which absorbs all scattered flelds. The 
mesh was constructed with the software built-in algorithm, which generates a free mesh 
consisting on tetrahedral elements. The subdomain representing the rod was constructed 
in all cases with mesh elements constrained to a maximum element edge size of 15 nm 
and a growth rate of 1.4, meaning that adjacent elements to a given one should not be 
bigger than 1.4 times the size of it. This allows to correctly map the near-fleld zone 
where the flelds decay rapidly. The rod is adaptively meshed with elements that are 
flner in the cups of the rod and coarser along the cylindrical body. As an illustrative 
example, in the case of a 340x30 nm rod, the subdomain representing the rod consisted 
on 18337 mesh elements, while the subdomains modeling the external medium and the 
PML consisted on 48783 mesh elements. The system was solved using PARDISO direct 
solver and involved, in this case, 427502 degrees of freedom. As the length of the rod 
was increased the number of mesh elements grew, and so did the number of degrees of 
freedom to be solved. The scattering efficiency was calculated as 



The integration was performed over a spherical boundary of radius Rcaic- We chose it to 
be the boundary between the embedding medium and the PML, thus Rcaic = Ri- Arod 
is the transverse geometrical area of the rod, which, in the case of normal incidence, is 
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simply 2RL + {it — 4) R"^ . Einc is the incident field amplitude, and Efar is the vectorial far- 
field of the scattered field obtained with the COMSOL implementation of the Stratton- 
Chu formula [36]. Please notice that, as previously pointed out by other authors [37] , 
COMSOL automatically introduces the factor R^^i^ when performing the integration, so 
that it has to be canceled in order to obtain the correct value for Qsca- 

For calculations of light scattering by an infinitely long nanowire in [6| we made use 
of our own implementation [381 EH] of Green's theorem surface integral equations, which 
are written for parametric surfaces describing particles with arbitrary shape by means 
of Giehs' "Superformula" [ID]. 

The dielectric function of Ag was modeled as a sum of Drude and Lorentz terms 
for all geometries. 



Appendix B. Details on fitting of scattering efficiency to Fano-like lineshape 

With respect to the slowly varying amplitude A{uj) in ([3]), we have heuristically chosen 



Aiuj) 



(B.l) 



where ujmin = 0.5wi, /? < 1 and Aq = |Ao|e*'^'. Best-fitting curves of Qsca to 
presented in figure |2| figure |3] and figure |5] were obtained by means of an iterative 
implementation of Levenberg-Marquardt algorithm [^ H3]. 

In order to place in context the information on the evolution of q parameter included 
in figure |2} we summarize in figure Bl all our fitting duties with respect to wi, ws, 61, 63 
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Figure B2. Obtained values of wi, W2, ^^3, ^3 for every angle a in figure [S] 



for all spheroids in figure [T] (L G [100,400] nm, D = 30 nm, Ed = 2.25). Logarithmic 
scale is used for the sake of readability. As expected from standard antenna theory, 
ci;i,a;3 scale linearly with 1/L. Interestingly, &1/63 ratio slightly decreases as L increases, 
thus fading the difference between "continuum-like" and "localized-like" resonances. 

For the case of p-polarized light impinging with oblique incidence on a single Ag 
spheroid surrounded by Ed = 2.25 (see figure |3]), parameters wi, a;2, 1^3, &2, ^3 remain 
almost invariant with respect to angle a, as can be seen in B2 This is consistent with 
the lack of interaction between resonances with different parity. 

Finally, we present in figure B3 the obtained ui,Us,bi, 63 for all spherocyhnders in 
figure [5] {L G [300,600] nm, D = 30 nm, Ed = 2.25). Except for numerical values, all 
parameters L follow the same trends than in Bl, thus confirming the validity of our 
Fano-like model when applied to a generic nanorod geometry. 
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